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A desirable feature of models with large extra dimensions and a TeV range string scale is the 
possibility to understand an ultralight sterile neutrino needed for a simultaneous understanding 
of solar, atmospheric and LSND results as a Kaluza-Klein state of a bulk neutrino. How these 
oscillation data are understood in quantitative detail is still not clear. We have recently suggested a 
new way to understand the solar neutrino data in this framework by a combination of vacuum and 
small angle MSW transition of v e to sterile KK modes of the bulk neutrino. This mechanism can 
be embedded into an economical extension of the standard model in the brane to understand the 
atmospheric and LSND data. After a brief discussion of how small neutrino masses arise naturally 
in extra dimensional models, we review this new suggestion. 
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03 ■ I. INTRODUCTION 
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One of the important predictions of string theories is the existence more than three space dimensions. For a long 
time, it was believed that these extra dimensions are small and are therefore practically inconsequential as far as low 
energy physics is concerned. However, recent progress in the understanding of the nonperturbative aspects of string 
| theories have opened up the possibility that some of these extra dimensions could be large [|]J|] without contradicting 
. observations. In particular, models, where some of the extra dimensions have sizes as large as a milli-meter and 
where the string scale is in the few TeV range have attracted a great deal of phenomenological attention in the past 
two years. The basic assumption of these models, largely inspired by generic observations in string theories is that 
the space time has a brane-bulk structure, where the brane is the familiar (3+1) dimensional space-time, with the 
standard model particles and forces residing in it and bulk consists of all space dimensions where gravity and other 
possible gauge singlet particles live. The main interest in these models is of course due to the fact that the low string 
scale provides an opportunity to test them using existing collider facilities. 

A major challenge to these theories comes from the neutrino sector, the first problem being how one understands 
O ' the small neutrino masses in a natural manner. The conventional seesaw || explanation which is believed to provide 
the most satisfactory way to understand this, requires that the new physics scale (or the scale of SU(2)r x U(1)b-l 
i-C ■ symmetry) be around 10 12 GeV or higher and clearly does not work for these models. So one must look for alternative 
, ways. The second problem is that if one considers only the standard model group in the brane, operators such as 
LHLH/M* could be induced by string theory in the low energy effective Lagrangian. For TeV scale strings this would 
r> | obviously lead to unacceptable neutrino masses. 

One mechanism suggested in Ref. Q] is to postulate the existence of one or more gauge singlet neutrinos, vb in the 
bulk which couple to the lepton doublets in the brane. After electroweak symmetry breaking this coupling can lead to 
neutrino Dirac masses, which are suppressed by the ratio M*/Mp£. This is sufficient to explain small neutrino masses 
and owes its origin to the large bulk volume that suppresses the effective Yukawa couplings of the Kaluza-Klein (KK) 
models of the bulk neutrino to the brane fields. In this class of models, naturalness of small neutrino mass requires 
that one must assume the existence of a global B-L symmetry in the theory, since that will exclude the undesirable 
higher dimensional operators from the theory. 

An alternative possibility [pi is to consider the brane theory to have an extended gauge symmetry which contains 
B-L symmetry as a subgroup, which will eliminate the higher dimensional operators. This is perhaps more in the 
spirit of string theories, which generally do not allow global symmetries. Phenomenological considerations, however, 
require that the local B — L scale and hence the string scale be of order of 10 10 GeV or so. The extra dimensions in 
these models could of course be large. 
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Regardless of which path one chooses for understanding small neutrino masses, a very desirable feature of these 
models is that if the size of extra dimensions is of order of a milli-meter, the KK excitations of the bulk neutrino have 
masses in the range of 10~ 3 eV, which is in the range needed for a unified understanding of oscillation data 0. 

In this paper, we will focus on TeV scale models and discuss their implications. In this class of models, one gets 
for the Dirac masses of the familiar neutrinos Ve.^.T 

hv wk M* 

The righthanded neutrino in the Dirac mass is a mode of the bulk neutrino i/g. For M* ~ 10 TeV, this leads to 
m v ~ 10 -4 /i eV. It is encouraging that this number is in the right range to be of interest in the discussion of solar 
neutrino oscillation if the Yukawa coupling h is appropriately chosen. 

If one takes these models for the neutrinos seriously, some immediate questions arise: first, while the KK modes 
of the bulk neutrinos are candidates for sterile neutrinos if the bulk radius is in the range of millimeters, are their 
small masses natural? Secondly, how well do the models of this type describe the observed oscillation data from solar, 
atmospheric and LSND experiments ? It is clear that they provide a tempting possibility to simultaneously explain 
all data since the needed sterile neutrino could be one of the KK modes of the bulk neutrino, as already emphasized. 

We discuss the answer to the first question in section 2. As far as the second question is concerned, several recent 
papers have addressed this issue ||j|-[ll[]. I n particular, in Ref. || it has been shown that while the overall features of 
the solar and atmospheric data can be accomodated in minimal versions of these models with three bulk neutrinos, it 
is not possible to simultaneously explain the LSND observation for the — v e oscillation probability, and one must 
incorporate new physics in the brane. 

On a phenomenological level, a first glance at the values of the parameters of the model such as mo from Eq.(l) 
and R^ 1 ~ 1CP 3 eV suggests that perhaps one should seek a solution of the solar neutrino data in these models 
using the small angle MSW mechanism |(|. However, present Super-Kamiokande recoil energy distribution seems to 
disfavor such an interpretation, although any definitive conclusion should perhaps wait till more data accumulates 
p2| . In any case if the present trend of the data near the higher energy region of the solar neutrino spectrum from 
Super-Kamiokande persists, it is likely to disfavor the small angle MSW solution and favor a vacuum oscillation. 
We therefore pursue the possibility of vacuum oscillation between the v e and v s to solve the solar neutrino puzzle 
motivated by a model that leads to desired parameters in a natural manner |13|j . 

As is known from many discussions in literature, chlorine results play a pivotal role in deciding on the nature of 
the oscillation solution to the solar neutrino puzzle. To fit water and the Gallium |]l9| data in conjunction with the 
Chlorine data, one must "kill" the 7 Be neutrinos almost completely. In the vacuum solution, it is in general hard to 
achieve this without simultaneously suppressing the pp neutrinos. Therefore, the general tactic adopted is to suppress 
the 7 Be neutrinos as much as possible (but not quite to zero) and suppressing the 8 B neutrinos to less than 30% 
or so. This, however, suppresses the higher energy 8 B contributions below the level observed by Super-Kamiokande 
collaboration. In the case of oscillation to active neutrinos, this can be fixed by including the neutral current scattering 
which contributes about 16% of the charged current data. On the other hand, this disfavors the vacuum oscillation to 
sterile neutrinos which do not have a neutral current component. As we will show later |l3|] , if the sterile neutrino is 
a bulk mode, its KK modes play a significant role in resolving this problem due to the existence of MSW transitions 
of v„ to these modes for values of extra dimension sizes of interest 



II. NATURALNESS OF ULTRA LIGHT STERILE NEUTRINOS IN BRANE-BULK MODELS 

Before getting to the detailed discussion of the model, we make some general comments about the naturalness of 
the ultralightness of the bulk neutrino. The bulk neutrino "self mass" terms are constrained by the geometry of the 
bulk and could therefore under certain circumstances be zero. If that happens, the only mass of the KK states of the 
vb will arise from the kinetic energy terms such as vb^ 1 ' di^s, where 1 = 5, 6... and will be given by n/R where R is 
the radius of the extra dimensions. In such a situation, an ultralight vb.kk arises naturally. 

The key to naturalness of the ultralight bulk neutrino is the geometry that forbids both Dirac and Majorana mass 
terms. Let us give a few examples. In five dimensions, if we impose the Zi orbifold symmetry (y — > —y), then it 
follows that the Dirac mass vanishes. Now if we impose lepton number symmetry in the brane, the Majorana mass 
vanishes, leaving us with no mass term for the bulk neutrino in 5-dimensions. 

Another interesting example is the 10-dimensional bulk, where the bulk neutrino is a 16-component spinor, which 
when reduced to 4-dimensions leads to eight 2-component spinors. The interesting point is that for a 16-dimensional 
spinor, one cannot write a Dirac or Majorana mass term consistent with 10-dimcnsional Lorentz invariancc. In this 
case, there is no need for assuming lepton number to get an ultralight sterile neutrino. A similar situation is also 
expected in six dimensions if we choose the bulk neutrino to be a 4-component complex chiral spinor. 
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III. BASIC IDEAS OF THE MODEL AND THE MIXING OF BULK MODES 



The model we will present will consist of the standard model in the brane (suitably extended to include some 
additional Higgs fields) and one bulk neutrino. The bulk may be five, six or higher dimensional; we will assume that 
only one of those extra dimensions is large. The bulk neutrino will be assumed to couple to vj^. Let us begin our 
discussion by focusing on the v e in the brane and one bulk neutrino. Obviously, the fields that could propagate in the 
extra dimensions are chosen to be gauge singlets. Let us denote bulk neutrino by v B (x tJl , y). It has a five dimensional 
kinetic energy term and a coupling to the brane field L(x fi ) given by 



C = KLHi> BR (x,y = 0) + dy v B L(x,y)d 5 v BR (x,y) + h.c, 



(2) 



where from the five dimensional kinetic energy, we have only kept the 5th component that contributes to the mass 
terms of the KK modes in the brane; H denotes the Higgs doublet, and k = h-M— the suppressed Yukawa coupling. It 
is worth pointing out that this suppression is independent of the number and radius hierarchy of the extra dimensions, 
provided that our bulk neutrino propagates in the whole bulk. For simplicity, we will assume that there is only one extra 
dimension with radius of compactification as large as a millimiter, and the rest with much smaller compactification 
radii. The smaller dimensions will only contribute to the relationship between the Planck and the string scale, but 
their KK excitations will be very heavy and decouple from neutrino spectrum. Thus, all the analysis could be done 
as in five dimensions. 

A second point we wish to note is that we will include new physics in the brane that will generate a Majorana mass 
matrix for the three standard model neutrinos as follows: 
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The origin of this pattern of brane neutrino masses will be discussed in a latter section. In this section let's discuss 
the mixing pattern of the bulk neutrinos with the brane ones. For this we will assume that mo 3> <%; as a result the 
v^ r in effect do not affect the mixing between the bulk neutrino modes and the v e . 

To discuss the effect of the bulk modes, note that the first term in Eq. (Q) will be responsible for the neutrino mass 
once the Higgs field develops its vacuum expectation. The induced Dirac mass parameter will be given by m = nv, 
which for M* = 1 TeV is about h ■ 10 -5 eV. Obviously this value depends only linearly on the fundamental scale. 
Larger values for M* will increase m proportionally. After introducing the expansion of the bulk field in terms of the 
KK modes, the Dirac mass terms in (|) could be written as 
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where our notation is as follows: v' B represents the KK excitations; the off diagonal term \f2m is actually an infinite 
row vector of the form y/2m(l, 1, • ■ •)• The operator 85 stands for the diagonal and infinite KK mass matrix whose 
rt-th entry is given by n/R. This notation was introduced in |^] to represent the infinite mass matrix in a compact 
manner. 

Using this short hand notation makes it easier to calculate the exact eigenvalues and the eigenstates of this mass 
matrix H . Simple algebra yields the characteristic equation 
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(5) 



where m n is the mass eigenvalue ^|||. The equation for eigenstates is 



1 In general, the bulk neutrino will couple to all three species of brane neutrinos; however, when the ^ MjT masses are much 
bigger than the mass parameters in the v E and vb sector, one can decouple the i^ jT and restrict oneself to only the Ve,VB 
sector, to analyse the effect of the bulk neutrino, as we do in the paper. 
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where the sum over the KK modes in the last term is implicit, and N n is the normalization factor given by 

2 (m„ - See) 2 



N 2 = l + m 2 ir 2 R z 



Using the expression (p) , we can write down the weak eigenstate vl in terms of the massive modes as 
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Thus, the weak eigenstate is actually a coherent superposition of an infinite number of massive modes. Therefore, even 
for this single flavour case, the time evolution of the mass-eigenstates involves in principle all mass eigenstates and is 
very different from the simple oscillatory behaviour familiar from the conventional two or three neutrino case. The 
survival probability depends strongly on the parameter £, reflecting the universal coupling of all the KK components 
of vb with vl in We will be interested in the small £ region of the parameter space. In this case it is easy to see 
that the coupling of the nth mode is given by £/n. 

Note that in the limit of S ee — 0, the v e and vb,o are two, two-component spinors that form a Dirac fermion with 
mass m. Once we include the effect of bee ^ 0, they become Majorana fermions with masses given by: m\ ps — (5 ee /2+m 



and 77i2 ~ — <5ee/2 — to, and they are maximally mixed; i.e., the two mass eigenstates are V\$ 
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Thus as the , 



produced in a weak interaction process evolves, it oscillates to the pb,o state with an oscillation length L ~ E/(2mS ee )- 
Later on in this paper we will find that for natural values of to, S ee in this model, L is of order of the Sun-Earth 
distance so that our model leads to vacuum oscillation of the solar neutrinos. 

Furthermore since the v e also mixes with the KK modes of the bulk neutrinos with a mass difference square of 
order 10 -6 eV 2 , there is MSW resonance transition of v e to vb.kk modes. In our discussion of the solar neutrino 
puzzle, we will take into account both these effects. 

In our model, we will also have too in the brane neutrino mass matrix which is of order eV so that we can account 
for LSND results. Atmospheric neutrino oscillation is now purely Va-v r oscillation. 



IV. SOLAR NEUTRINO DATA BY A COMBINATION OF VACUUM AND MSW OSCILLATIONS 

Let us now discuss how the solar neutrino data is explained within this model. To discuss the MSW effect for 
the case of bulk neutrinos, we need to include the matter effect in the diagonalization of the infinite dimensional 
neutrino mass matrix. The eigenvectors and eigenvalues of this matrix can be found in the presence of matter, when 
the squared mass matrix, M 2 , is replaced with M 2 + 2EH\, where H\ is p e — \/2GF{n e — 0.5n n ) when acting on v e , 
and is zero on sterile neutrinos. Define 

E P e | I | 1 E P e \\ 

m k See V \m k S ee J 
w k = 1 in vacuum. The characteristic equation becomes 

Tmf_ TTTOfc 

777fc = W k dee H Cot , (10) 
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The eigenvectors are as in Eq. except the coefficients of vqb and v' B _ acquire an additional factor l/w n , and the 

normalization becomes N 2 = 1 H I n 2 m 2 R 2 + ™% ) ™ n ■ 

In fitting the solar data the simplest way to reconcile the rates for the three classes of experiments is to "kill" the 
7 Be neutrinos, reduce the 8 B neutrinos by half and leave the pp neutrinos alone. To achieve this in the VO case, one 
may put a node of the survival probability function P ee around 0.86 MeV. However, for an arbitrary node number, the 
oscillatory behaviour of P ee before and after 0.86 MeV cannot in general satisfy the other two requirements mentioned 
above. If one uses the first node to "kill 7 Be" , then for 8 B neutrino energies the P ee is close to one and not half 
as would be desirable; on the other hand, if one uses one of the higher nodes (higher Am 2 ), then pp neutrinos get 
reduced. The strategy generally employed is not to "shoot" for a node at the precise 7 Be energy but rather somewhat 
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away so that it reduces 7 Be to a value above zero. This requires that one must reduce the 8 B neutrinos by much 
more than 50%, so one can fit Chlorine data. The water data then requires an additional contribution, which, in the 
case of active VO, is provided by the neutral current cross section amounting to about 16% of the charged current 
one. Thus in a pure two-neutrino oscillation picture, VO works for oscillation to active neutrinos but does not work 
for active to sterile oscillation. It is here that the large extra dimensions come to the rescue as we see below, thereby 
keeping the active to sterile VO fit to solar neutrino data viable. 

To see this note that in our model both vacuum oscillations and MSW oscillations are important. This is because 
the lowest mass pair of neutrinos is split by a very small mass difference, whereas the KK states have to be separated 
by < 10 -2 eV because of the limits from gravity experiments. We can then use the first node of P ee to suppress the 
7 Be. Going up in energy toward 8 B neutrinos, the survival probability, which in the VO case would have risen to 
very near one, is suppressed by the small-angle MSW transitions to the different KK excitations of the bulk neutrino. 
This is the essence of our new way to fit the solar neutrino observations [^3| . 

To carry out the fit, we studied the time evolution of the v e state using a program that evolved from one supplied 
by W. Haxton jlj] . The program was updated to use the solar model of BP98 [^5) and modified to do all neutrino 
transport within the sun numerically. For example, no adiabatic approximation was used. Changes were also necessary 
for oscillations into sterile neutrinos and to generalize beyond the two-neutrino model. Up to 16 neutrinos were allowed, 
but no more than 14 contribute for the solutions we considered. 
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FIG. 1. Energy dependence of the v e survival probability when R — 58jum, mR = 0.0093, 8 ee — 0.84 x 10 -7 eV. The 
dot-dashed part of the curve assumes the radial dependence in the Sun for neutrinos from the pp reaction, the solid part 
assumes ls O radial dependence, and the dashed part assumes S B radial dependence. 

For comparison with experimental results, tables of detector sensitivity for the Chlorine and Gallium experiments 
were taken from Bahcall's web site The Super-Kamiokande detector sensitivity was modeled using |Iq |, where 
the percent resolution in the signal from Cerenkov light, averaged over the detector for various total electron energies, 
is provided. For details see 

Calculations of electron neutrino survival probability, averaged over the response of detectors, were compared 
with measurements. While theoretical uncertainties in the solar model and detector response were included in the 
computation of \ 2 as described in Ref. p7| , the measurement results given here include only experimental statistical 
and systematic errors added in quadrature. The Chlorine survival probability, from Homestake [jUst , is 0.332 ± 0.030. 
Gallium results |l| for SAGE, GALLEX and GNO were combined to give a survival probability of 0.579 ± 0.039. 
The 5.5 — 20 MeV 1117 day Super-K experimental survival probability [jl9| is 0.465 ± 0.015. The best fits were with 
R w 58/iTO, mR around 0.0093, and 5 ee ~ 0.84 x 10 -7 eV, corresponding to 6m 2 ~ 0.53 x 10 -11 eV 2 . These parameters 
give average survival probabilities for Chlorine, Gallium, and water of 0.386, 0.533, and 0.460, respectively. They give 
a v e survival probability whose energy dependence is shown in Fig. |l| For two-neutrino oscillations, the coupling 
between v e and the higher mass neutrino eigenstate is given by sin 2 28, whereas here the coupling between v e and the 
first KK excitation replaces sin 2 29 by 4m 2 R 2 — 0.00035. 

Vacuum oscillations between the lowest two mass eigenstates nearly eliminate electron neutrinos with energies of 
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.63MeV/(2n+l) for n = 0, 1, 2, ... . Thus Fig. [TJ shows nearly zero v e survival near 0.63 MeV, partly eliminating 
the Be contribution at 0.862 MeV, and giving a dip at the lowest neutrino energy. Increasing S ee moves the low 
energy dip to the right into Gallium's most sensitive pp energy range, making the fit worse. Decreasing S ee increases 
Gallium, but hurts the Chlorine fit by moving the higher energy vacuum oscillation dip further to the left of the 7 Be 
peak. MSW resonances start causing the third and fourth eigenstates to be significantly occupied above ~ 0.8 MeV, 
the fifth and sixth eigenstates above ~ 3.7 MeV, the 7'th and 8'th above ~ 8.6 MeV, and the 9'th and 10'th above 
~ 15.2 MeV. Fig. [I] shows dips in survival probability just above these energy thresholds. The typical values of the 
survival probability within the a B region (~ 6 to ~ 14 MeV) are quite sensitive to the value of mR. As can be seen 
from Eq. [?], higher mR increases 1/N « m/m n « mR/n for various n, and thereby increases v e coupling to higher 
mass eigenstates, strengthens MSW resonances, and lowers v e survival probability. 

The expected energy dependence of the v e survival probability is compared with preliminary Super-K data |[9[ in 
Fig. ^. The uncertainties are statistical only. The parameters used in making Fig. || were chosen to provide a good 
fit to the total rates only; they were not adjusted to fit this spectrum. But combining spectrum data with rates using 
the method described in Ref. p0| , with numbers supplied by Super-K |2J, gives \ 2 — 15.6 ("probability" 74%) for 
the spectrum predicted from the fit to total rates. It's incorrect to calculate probability as if there were (number of 
data points) - (number of parameters) = 17 degrees of freedom, but if we do, x 2 corresponds to a "probability" of 
55%. 

One may also seek fits with 8 ee constrained to be very small, thereby eliminating vacuum oscillations. The best such 
fit had x 2 — 5-5 ("probability" 14%). The same parameters then used with the Super-K spectrum gave \ 2 — 18.7 
("probability" 35%). 
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FIG. 2. Super-Kamiokande energy spectrum: measured [19] preliminary results based on 1117 days (error bars) and predicted 
(curve) for the same parameters as in Fig. 1. The curve is not a fit to these data. 

The seasonal effect was computed for a few points on the earth's orbit. If r is the distance between the earth and 
the sun, 

— = l + ecos(8-e ), (11) 
r 

where r is one astronomical unit, e = 0.0167 is the orbital eccentricity, and 9 — 9q rs 2ir(t — f ), with t in years and 
<o = January 2, 4h 52m. Table ^ shows very small seasonal variation. 

To understand the mass pattern used, consider L e +Lfj,—L T symmetry for neutrinos, with L e = 1 for i>b- The allowed 
mass pattern for (y e , v^, v T ) is given by Eq. |^ with all J's (except S eT ^ T ) set to zero. The remaining S entries arise after 
we turn on the symmetry breaking. For an explicit realization, we augment the standard model by the singlet charged 
Higgs rj + , h ++ which are blind with respect to lepton number and SU(2)l triplet fields, A eifJl . T , with Y = 2 which carry 
two negative units of lepton numbers L ei/J/iT , respectively. The Lagrangian involving these fields consists of two parts: 
one Cq which is invariant under (L e + — L T ) number and contains terms (rjL^L T ,riL e L T ^ h ++ e~^r^, h ++ fi~^r^) 
and L e L e A e , L^L^A^, L T L T A T and a soft breaking term £ a = h ++ '(£)j =ej|ijT M U A? + Af pA p A T + Af 0r A e A T ) + h.c. 
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With these couplings, the neutrino Majorana masses arise from two-loop effects (similar to the mechanism of ref. 
p2|), and we have the L e + — L T violating entries &y ~ c m ei m ej . Using S ee ~ 10~ 8 eV, then 5 TT ~ 1CP 2 , as 
would be required to understand the atmospheric neutrino data. 

Finally, models with bulk sterile neutrinos lead to new contributions to low energy weak processes, which have been 



addressed in several papers |23|. In the domain of astrophysics, they lead to new contributions to supernova energy loss, 
as well as to the energy density in the early universe, which can influence the evolution of the universe. Currently these 
issues are under discussion pi|Jq^q], and if neutrino oscillation data favor these models, any cosmological constraints 
must be addressed. Note, however, that this model is less sensitive to these issues than other fits because for VO 
Am 2 is an order of magnitude smaller, for MSW sin 2 29 is more than an order of magnitude smaller, and there is 
only a single KK tower starting from a very small mass |24j] . 

TABLE I. Predicted seasonal variations in u e fluxes, excluding the 1/r 2 variation. The model assumed fio = 0.32 x 10 -2 eV, 
m = 0.34 x 10~ 4 eV, and S ee = 0.78 x 10 -7 eV. 
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